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We study the frustrated dimer-plaquette quantum spin chain for ferromagnetic 
dimer bonds. This quantum system undergoes a series of first-order ground-state 
phase transitions driven by frustration or by a magnetic field. We find that the 
different nature of the ground-state phases has a strong influence on the magne- 
tization curve as well as the low-temperature thermodynamics. In particular, the 
magnetization curve exhibits plateaus and jumps, the number of which depends on 
the strength of frustration. The temperature dependence of the susceptibility may 
either show activated spin-gap behaviour for small frustration or Curie like param- 
agnetic behaviour for strong frustration. 

PACS numbers: 75.10.Jm 

INTRODUCTION 

Over the last years much attention has been concentrated on the physics of low- 
dimensional quantum spin systems. In particular, zero-temperature phase transitions and 
exotic magnetization curves in frustrated quantum magnets are in the focus of investiga- 
tions. Besides continuous quantum phase-transitions also remarkable first-order transitions 
can be driven by frustration. The theoretical studies have benefitted from recent experimen- 
tal results on low-dimensional spin-half antiferromagnets like quasi-onedimensional ladder, 
zigzag and spin-Peierls systems as well as gapped quasi-t wo dimensional quantum magnets 
like CaV.O, i i and SrCu2{B02h i i §• 

Among the various models for low-dimensional quantum magnets the spin-half frustrated 
dimer-plaquette chain (FDPC) has been discussed in a series of papers in recent years. This 
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FIG. 1: The fustrated dimer-plaquette chain (FDPC). 
FDPC was introduced in and its Hamiltonian reads (see also Fig. |I]) 

Np Np Np 

H = H,, + Hf = J,Y. S^S^ + E(S: + SD(S^ + Sr^) + J/ E s^s^ (1) 

n=l n=l n=l 

where A'p = N/ 4 is the number of plaquettes (unit cells) and is the number of spins. We 
consider finite chains of spins with periodic boundary conditions. The Hamiltonian fulfills 
the important relation [H, (S^ + S^)^]_ = 0. Hence the vertical dimer spins S" and form 
a composite spin = (S^J + S;^) with eigenvalues {S]^^,f = 5^^(55, + 1) and = {1,0}. 
The physics of this model was discussed in the literature for antiferromagnetic couplings 

In the unfrustrated version {Jf = 0, H = Hdp) the model may serve as the onedimensional 
counterpart of the 1/5-depleted square lattice of CaV^Og |rT|. In this case (J/ = 0) 

all composite spins have eigenvalues S*"^ = 1 in the singlet ground state as well as in the 
first triplet excitation and the ground-state physics of (|^) is determined by the competition 
of dimer (J^) and plaquette (Jp) bonds. In the limit ^ Jp the low-energy physics of 
the model is that of the Haldane chain with an effective coupling J^jj = Jp/2Jd, whereas 
for Jp Jd a plaquette phase is realized 0. The singlet-triplet excitation gap is always 
finite for Jp > 0. A special feature of the model is the existence of a class of exact product 
eigenstates, for which the composite spins of certain plaquettes form a vertical dimer singlet 
(i.e. S^fj = for certain n). The finite strips (fragments) between two vertical dimer singlets 
contain vertical dimer triplets S*"^ = 1 and are decoupled from each other. Therefore these 
product eigenstates correspond to a fragmentation of the chain. 

Taking into account frustration Jf > the formation of vertical dimer singlets becomes 
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energetically more favourable and the system undergoes at a critical frustration Jj — 
f{Jd, Jp) a first-order quantum phase transition from the collective ground state with S^^ = 1 
for all n to a fragmented dimer product ground state 

i^o...o) = n 2-^ (I - 1 CK)) n Kp), (2) 

n=l n=l 

with 5""^ = for all n. For > the dimer state \d^p) ) is a singlet and the energy of 
|^o...o) is E/Np = -SJf/A - ^Jd/^. 

In the special limit Jj = Jd the critical value is Jj"* = 1. 2210 Jp [|]. In this limit the FDPC 
has a close relation to the two-dimensional spin model for SrCu2 (603)2 @, and is called 
orthogonal-dimer chain. Similar to the quasi-t wo dimensional spin system SrCu2 (603)2 P, |^ 
the FDPC exhibits nontrivial magnetization plateaus [§, P]. A novel property of the FDPC, 
not found so far in quantum spin systems, is the existence of an infinite sequence of plateaus 

i- 

Very recently, the FDPC has been considered for arbitrary spin quantum numbers s > 
1/2 and a series of 2s first-order quantum phase transitions has been found by Koga and 
Kawakami [^. Koga and coworkers argue that the first-order phase transitions described 



by the FDPC possess most of the essential features inherent in frustrated quantum spin 
systems [P, 10|. In that sense the FDPC may serve as a prototype model for spin sytems 



showing first-order quantum phase transitions. 

In the present paper we extend the model to ferromagnetic bonds Jd- We notice, that 
the fragmentation of the chain due to vertical dimer singlets occurs for Jd < 0, too. In 
particular, the dimer product state @ is an eigenstate of (|l|) also for J^ < 0, however, the 
dimer state Id^^j) in |\I'o...o) is a triplet. The energy of |^o...o) is then E /Np = —3Jf /4: + Jd/4:. 

GROUND STATE PHASE DIAGRAM 

We start from the energy eigenvalues E of model (|l|) and write the dependence of E on 
Jf in explicit form 

E{Jp, Jd, Jf) = EdpiJp, Jd) + Jf (^iV, - ^Ns^ , (3) 

where Ns is the number of vertical dimers with composite spin 5*"^ = and Nt is the number 
of vertical dimers with composite spin S'^b = 1- We have Ns + Nt = Np. Since every vertical 
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FIG. 2: Change of the ground state with increasing frustration Jf. For small fustration Jf all 
composite spins are in the triplet state and there is no fragmentation. Increasing Jj some composite 
spins may prefer a vertical singlet state (indicated by vertical lines) and the chain may split into 
identical fragments of length k (the spin-spin correlation is zero along the Jp bonds neighbouring 
a vertical singlet as indicated by dotted lines). For large Jf all composite spins are in the singlet 
state and the dimer product state (0) is the ground state. 



dimer singlet leads to an energy gain of Jf in the second part of the energy the frustration 
favours states with S^^, = 0, i.e. a fragmentation of the chain. The lowest state of these 
finite fragments has a total spin Sfrag = 1 for J^ < 0. We define the length A; of a fragment, 
as the number of triplet composite spins 5*"^, = 1 between two singlet composite spins. 
According to the so-called linear programming scheme the extrema of the energy belong 



to states with fragments of identical length k . In what follows we use the notation <k> 
for a fragmented state consisting of identical fragments of length k. The dimer product 
state (0) corresponds to fragments of length k = 0, i.e. |\E'o...o) =<0>. It is evident that the 
system undergoes at least one transition from the nonfragmented ground state <oo> with 
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5""^ = 1 (n = 1 . . . Np) at small Jj to the dimer product ground state <0> with S*"^ = 
{n = 1 . . . Np) for large Jf (see Fig. 0) . As mentioned above, for antiferromagnetic Jd indeed 
we have a direct transition <oo>^<0>. However, for ferromagnetic Jd the question for 
the existence of intermediate ground states <k> with fragments of finite length < A; < cx3 
needs more attention. 

We start from the ground state <oo> for low frustration. The energy per unit cell of this 
state is 

= EUNp = + J//4, (4) 

where Coo is the energy of the ground state for the unfrustrated chain {Jf = 0). We determine 
Coo by finite-size extrapolation of chains of length N = 8, 16, 24, 32 as well as perturbation 
theory. The energy per unit cell e{ of a fragmented state <fc> is given by the sum of the 
energies + kJf/i of the fragments and of the energies SJf/A of the singlets separating 
the fragments 

f _ E, + {k~ 3)Jf/4 

Ek is the energy of the unfrustrated ( J/ = 0) fragment of length k consisting of k plaquettes 
and k + 1 dimer bonds Jd and corresponds to the energy of a finite dimer-plaquette chain 
with 4/c + 2 spins and open boundary conditions. 

Using Lanczos algorithm we can exactly caculate the energy of the chain fragments E^ 
up to = 8 (i.e. = 34). For larger fragment lengths and for the state < oo > we 
need approximations. In the limits |t/d|/|Jp| <^ 1 and |t/p|/|t/d| <^ 1 perturbation theory is 
appropriate since for Jd = and Jp = the ground states are known as simple product states. 
We use here Rayleigh-Schrodinger perturbation theory up to 3rd order. The comparison of 
perturbation theory with Lanczos results shows that the perturbation theory yields reliable 
results even for \Jd/ Jp\ ~ 1- Fig. |] shows the lowest energies of the state <k> with 
/c = 0, 1, 2, 3, . . . und k = oo versus Jf for several Jd- The intersection points of the lowest 
lines determine the positions of first-order quantum phase transitions. It becomes obvious 
that the number of phase transitions changes from one to four increasing the strength of 
ferromagnetic Jd- The existence of further transitions can be excluded by analyzing the 
^-dependence of E^ according to Niggemann et al. [jl2|. The transition point Jj^'^^ between 
two ground-state phases <ki> and <k2> are obtained by the relation = e{2- Using Eqs. 
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FIG. 3: Energies e{ (Eq. (||)) versus Jj for various strengths of ferromagnetic dimer bonds 
{Jd = —0.1, Jd = —6, Jd = —100 and Jd = — oo) and fixed Jp = 1. For better comparison we 
have subtracted e{Jp = 0, J/ = 0) = Jd/4 from e{. Curves are presented for A; = 0, . . . , 7 (exact 
results) and for k = oo (finite-size extrapolation). The accuracy of the k = oo results corresponds 
to the thickness of the solid line. We have labeled by <A;> only the most relevant curves with 
k = 0,1, 2, 3, oo. The inset in the figure for Jd = —oo shows the transition region between <oo>, 
<3> and <2> with an enlarged scale. 

(I) and (I) we find 

j'/°^ = E,-{k + l)eoo ; j'}-''' = {k + l)Ek-i~{k)Ek. (6) 

The corresponding phase diagram is shown in Fig. ^. For completeness and comparison we 
have reconsidered the case of antiferromagnetic Jd P, |^ including perturbation theory and 
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FIG. 4: Ground-state phase diagram for ferromagnetic Jd < and fixed Jp = I (see text), 
solid line - perturbation theory, crosses - exact diagonalization 




FIG. 5: Ground-state phase diagram for antiferromagnetic Jd > and fixed Jp = I (see text), 
solid line - perturbation theory, crosses - exact diagonalization 
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an enlarged exact-diagonalization data set, see Fig. ^. For small strength of dimer coupling 
\Jd\ ^3 the situation is similar for both ferro- and antiferromagnetic J^, i.e. we have a 
direct transition <oo>^<0>. Increasing the strength of ferromagnetic dimer coupling to 
Jd ~ —4 the situation is changed and we find an additional intermediate phase <1>. Further 
increasing \ Jd\ we find another intermediate phase <2>, and finally we have a sequence of 
transitions < oo >^< 3 >— 2 >— 1 >— >< > at very strong ferromagnetic Jd- In this 
extreme limit ^ — oo the FDPC can be mapped onto the frustrated diamond chain and 



our results are in agreement with the results reported in Ref. |12 



INFLUENCE OF A MAGNETIC FIELD ON THE GROUND-STATE PHASES 

As discussed in Refs. ^ ^ the FDPC with antiferromagnetic Jd exhibits nontrivial mag- 
netization plateaus. Therefore we consider now the Hamiltonian (|l|) including an external 
magnetic field 

H = Hdp + Hf - hS^^t, (7) 

where the z-component of the total spin is given by S^^^ = J2i Si - For ferromagnetic Jd we 
have a more complex zero-field ground-state phase diagram and we can expect interesting 
effects caused by a magnetic field. One important difference to the antiferromagnetic case 
Jd > consists in the circumstance, that for ferromagnetic Jd any fragment of length k 
carries a finite total spin Sfrag = 1 in its lowest state for h = 0, whereas the total spin Sfrag 
for Jd > is zero for even k. Hence the fragmentation is favoured by the magnetic field, 
whereby short fragments are more favourable than long fragments. As discussed already 
above, the extrema of the energy belong to states with fragments of identical length k- For 
Jd < in zero field the fragmented state <k> {k finite) can carry total spin S^^f- = N/ {Ak+A). 
From the numerical inspection of the excitation spectrum of finite fragments we argue that 
their excited states with Sfrag > 1 are well separated from the lowest state with Sfrag = 1- 
Hence we can expect for parameter points in the zero-field phase diagram (Fig. ^) not too 
far from a transition line that the first-order transitions can also be driven by an external 
magnetic field. 

To be more precise we have studied the ground-state phases of a FDPC in a magnetic 
field with = 32 sites in the limit Jd —>■ — oo. The corresponding phase diagram is shown in 
Fig. ^. We mention, that for this chain length the fragmented state <2> does not fit to the 



9 




0.5 1 1.5 2 



Jf 



FIG. 6: Ground-state phase diagram for the FDPC of = 32 spins in a finite magnetic field h 
with Jp = 1 and large ferromagnetic Jd — > — oo. 

periodic boundary conditions and is therefore missing. (A finite-size calculation including 
all zero-field ground-state phases would require N = 48 sites which is currently beyond the 
available computer facilities.) We argue that except for missing the <2> phase the other 
ground-state phases should be well described by the = 32 system, since the correlation 
length in all phases is rather small. In particular, the transition lines between m = 1/4 
and m = 1/2 and between m = 1/2 and m = 1 presented in Fig. ^ are exact, since these 
ground-state phases belong to product states correctly described for = 32. 

Similar to the zero-field phase diagram there is no fragmentation for small frustration. 
However, supported by the field the fragmention sets in already at Jf ~ O.TlJp for h ^ 0.13 
instead of Jf ~ 1.30Jp for h = 0. In general, the phase transitions shown in Fig. ^ are 
shifted to lower values of frustration Jf. 

The ground-state phase diagram in finite magnetic field leads to interesting magnetization 
curves (Fig. |^). We define the magnetization as m = 2Sf^^/N, m is zero in a singlet state but 
unity in the fully polarized ferromagnetic state. Then the fragmented ground-state phases 
<oo>, <3>, <2>, <1>, <0> in Fig. ^ correspond to m = 0, m = 1/8, m = 1/6, m = 1/4 
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FIG. 7: Magnetization curves m{h) for the FDPC of = 32 spins with Jp = 1, large ferromagnetic 
Jd ^ —oo and various strengths of frustration Jf. 

and m = 1/2, respectively. 

For small Jf ^ 0.70 Jp we have a magnetization curve m{h) typical for unfrustrated 
gapped spin systems, like the Haldane chain [|^, |15|, |l^]. Due to the gap the curve starts 
with a plateau at zero field but then it goes continuously to saturation. For the finite system 
considered the continuous part is staircase like with small steps (see thin solid line in Fig. 
^ which is related to the corresponding line of a = 16 Haldane chain. For comparison 
we added the DMRG data for the Haldane chain of 60 sites (see thick solid line in Fig. ^ 
taken from [|15[. For intermediate frustration 0.71Jp ^ Jf ^ l.OSJp we have jumps, plateaus 



and continuous parts in the m{h) curve. However, for large enough Jf > l.OSJp we have 
no continuous parts in the curve, but only plateaus connected by jumps. For instance, for 
Jf = 1.2Jp (dotted line in Fig. 0) we have the same sequence of phase transitions as in Fig. 
^, but now driven by increasing the magnetic field from /i = 0to/i~0.05. As a result the 
magnetization m jumps from m = 0tom = l/8tom = l/6tom = l/4tom = 1/2 and 
further increasing h to saturation m = 1 (notice that m = 1/6 is missing for A^ = 32 in Fig. 
1^. For large frustration Jf < 1.5Jp we have an extreme m{h) curve consisting only of one 
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plateau at m = 1/2 followed by a jump to saturation m = 1. 

We remark, that the jump to the saturation m = 1 present for Jf ^ 0.97Jp is related to 



independent magnon excitations versus the fully polarized ferromagnetic state [0. 



We notice, that our result is in accordance with the general rule of Oshikawa et al. |T8 
that n{s — ffi) has to be an integer. In our case the period of a fragmented ground state 
<k> is Ak + 4, the spin s is one half and the magnetization per site fh is 1/(4A; + 4) (the 
magnetization per site ffi corresponds to one half of the magnetization m used so far). 



FINITE TEMPERATURES 



In this section we discuss some consequences of the ground-state phase diagram on the 
low-temperature thermodynamics. The different nature of the nonfragmented phase <cxd> 
and the fragmented phases <k> {k finite) has an important impact on thermodynamics. 
While the the ground state <oo> is a singlet with gapped triplet excitations, the fragmented 
ground states <k> consist of independent paramagnetic entities. As a consequence, in 
the former case both the susceptibility x the specific heat c are thermally activated 
and decay exponentially to zero if temperature T goes to zero. On the other hand, for a 
fragmented ground state the susceptibility shows a paramagnetic Curie like divergency for 
T — i> 0. Hence, varying Jf the low-temperature behaviour of x changes basically crossing the 
transition line between <cx3> and <3>. The behaviour for the specific heat near the phase 
boundaries can be more complex, since close to these boundaries we have quasi-degenerated 
low-lying levels which may lead to additional low-temperature peaks in c. 

To illustrate this behaviour we present in Fig. ^ and Fig. ^ the susceptibility and the 
specific heat for a FDPC of = 16 sites obtained by full diagonalization of the Hamiltonian 
matrix. The transition from spin-gap to paramagnetic behaviour is clearly seen in Fig. ^ 
between Jf = 1.2 and Jf = 1.4. In the same parameter region the specific heat shows a 
characteristic double-peak structure due to the quasi-degeneracy of states near the transition 
line. This effect should survive in the thermodynamic limit. 
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FIG. 8: Magnetic susceptibility % of the FDPC of A?^ = 16 spins with Jp = 1, ferromagnetic 
Jd = —400 and various strengths of frustration Jf. 



a 




FIG. 9: Specific heat c of the FDPC of iV = 16 spins with Jp = 1, ferromagnetic Jd = —400 and 
various strengths of frustration Jf. 
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SUMMARY 

In conclusion, we have found a series of first-order quantum phase transitions in a frus- 
trated one-dimensional quantum spin system which is closely related to the possibility of 
fragmentation of the considered chain. This series of phase transitions can be driven either 
by frustration or by a magnetic field. We emphasize that this fragmentation and the related 
consequences for physical properties are pure quantum effects not present in classical spin 
systems. The existence of different ground states has a strong impact on the magnetization 
curve and the low-temperature thermodynamics of the spin system. In dependence on the 
frustration the magnetization curve shows plateaus, jumps as well as continuous parts. 

Though the possibility of fragmentation has been observed also for frustrated two-leg 
spin ladder |15|, |l9l, in the spin ladder system only very simple ground states seem to be 
relevant for its magnetization curve. 

The low-temperature thermodynamics can be either spin-gap like for small frustration or 
paramagnetic for large frustration. Furthermore an additional low-temperature peak in the 
specific heat can appear. 

Acknowledgments 

This work was supported by the DFG (Ri 615/6-1). The authors are indebted to 
N.B. Ivanov for fruitful discussions and to A. Honecker for critical reading the manuscript 
and sending the DMRG data used in Fig. 0. 



[1] S. Taniguchi, T. Nishikawa, Y. Yasui, Y. Kobayashi, M. Sato, T. Nishioka, M. Kontani, and 

K. Sano, J. Phys. Soc. Jpn. 64, 2758 (1995). 
[2] T. Ohama, H. Yasuoka, M. Isobe, and Y. Ueda, J. Phys. Soc. Jpn. 66, 23 (1997). 
[3] H. Kageyama, K. Yoshimura, R. Stern, N.V. Mushnikov, K. Onizuka, M. Kato, K. Kosuge, 

CP. Slichter, T. Goto, Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 
[4] S. Miyahara and K. Ueda, Phys.Rev.Lett. 82, 3701 (1999). 

[5] K. Onizuka, H. Kageyama, Y.Narumi, K.Kindo, Y. Ueda, and T. Goto, J. Phys. Soc. Jpn. 
69, 1016 (2000). 



14 



[6] N.B. Ivanov, and J. Richter, Phys.Lett.A 232, 308 (1997). 

[7] J. Richter, N.B. Ivanov, and J. Schulenburg, J. Phys.: Condens. Matter 10, 3635 (1998). 

[8] A. Koga, K. Okunishi, and N. Kawakami, Phys. Rev. B 62, 5558 (2000). 

[9] J. Schulenburg and J. Richter, Phys. Rev. B 65 054420 (2002). 

[10] A. Koga and N.Kawakami, Phys. Rev. B 65, 214415 (2002). 

[11] N. Katoh and M. Imada, J. Phys. Soc. Jpn. 64, 4105 (1995). 

[12] H. Niggemann, G. Uimin, and J. Zittartz J. Phys.: Condens. Matter 9, 9031 (1997). 

[13] K. Ueda and S. Miyahara, J. Phys.: Condens. Matter 11, L175 (1999). 

[14] T. Sakai and M. Takahashi. Phys. Rev. B 43, 13383 (1991). 

[15] A. Honecker, F. Mila and M. Troyer, Eur. Phys. J. B 15, 227 (2000). 

[16] M. E. Zhitomirsky and I. A. Zahznyak, Phys. Rev. B 53, 3428 (1996). 

[17] J. Schulenburg, A. Honecker, J. Schnack, J. Richter and H.J. Schmidt, Phys. Rev. Lett. 88, 

167207 (2002). 

[18] M. Oshikawa, M. Yamanaka, I. Affleck, Phys. Rev. Letter 78, 1984 (1997). 

[19] F. Mila, Eur. Phys. J. B 6, 201 (1998). 



